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Abstract 



> 

m 

\ Krawtchouk matrices have as entries values of the Krawtchouk 

■ polynomials for nonnegative integer arguments. We show how they 

, arise as condensed Sylvester-Hadamard matrices via a binary shuffling 

\ function. The underlying symmetric tensor algebra is then presented. 

To advertise the breadth and depth of the field of Krawtchouk 
polynomials/matrices through connections with various parts of math- 
' ematics, some topics that are being developed into a Krawtchouk En- 

^ • cyclopedia are listed in the concluding section. Interested folks are 

encouraged to visit the website 



jhttp : //chanoir . math . siu . edu/wiki/KravchukEncyclopedia 



which is currently in a state of development. 

1 What are Krawtchouk matrices 

Of Sylvester-Hadamard matrices and Krawtchouk matrices, the latter are 
less familiar, hence we start with them. 

Definition 1.1 The N^^-order Krawtchouk matrix K^^^ is an (A^ + 1) x 
(A^ + 1) matrix, the entries of which are determined by the expansion: 

N 

{l + vf~Hl-vy = Y,v'KiP (1.1) 

i=0 
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Thus, the polynomial G{v) = (1 + v)^~^ (1 — vY is the generating function 
for the row entries of the j^^ column of K^^\ Expanding gives the explicit 
values of the matrix entries: 

-r^i:(-)'(0(r-7)- 

where matrix indices run from to N . 



Here are the Krawtchouk matrices of order zero, one, and two: 



= [ 1 ] 



1 1 
1 -1 



111 

2 0-2 
1 -1 1 



The reader is invited to see more examples in Table 1 of the Appendix. 



The columns of Krawtchouk matrices may be considered generalized bi- 
nomial coefficients. The rows define Krawtchouk polynomials: for fixed order 
iV, the i^^ Krawtchouk polynomial takes its corresponding values from the i^^ 
row: 

H3.N) = k\P (1.2) 

One can easily show that ki{j,N) can be given as a polynomial of degree 
i in the variable j. For fixed N, one has a system of + 1 polynomials 
orthogonal with respect to the symmetric binomial distribution. 

A fundamental fact is that the square of a Krawtchouk matrix is propor- 
tional to the identity matrix. 

(irW)2 = 2^-/ 

This property allows one to define a Fourier-like Krawtchouk transform on 
integer vectors. For more properties we refer the reader to [12]. In the 
present article, we focus on Krawtchouk matrices as they arise from corre- 
sponding Sylvester-Hadamard matrices. More structure is revealed through 
consideration of symmetric tensor algebra. 

Symmetric Krawtchouk matrices. When each column of a Krawtchouk 
matrix is multiplied by the corresponding binomial coefficient, the matrix 
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becomes symmetric. In other words, define the symmetric Krawtchouk 
matrix as 

where B^^^ denotes the (A^ + 1) x (A^ + 1) diagonal matrix with binomial 
coefficients, B^f ^ = (^), as its non-zero entries. 

Example. For = 3, we have 

13 3 1 " 

3 3-3-3 

3-3-3 3 

1-3 3-1 



Some symmetric Krawtchouk matrices are displayed in Table 2 of the 
Appendix. A study of the spectral properties of the symmetric Krawtchouk 
matrices was initiated in work with Fitzgerald |llj . 

Background note. Krawtchouk's polynomials were introduced by Mikhail 
Krawtchouk in the late 20's [TTJ [18]. The idea of setting them in a matrix 
form appeared in the 1985 work of N. Bose [2] on digital filtering in the 
context of the Cayley transform on the complex plane. For some further 
development of this idea, see [12] . 

The Krawtchouk polynomials play an important role in many areas of 
mathematics. Here are some examples: 

• Harmonic analysis. As orthogonal polynomials, they appear in the 
classic work by Szego [21]. They have been studied from the point of 
view of harmonic analysis and special functions, e.g., in work of Dunkl 
[SI [H] . Krawtchouk polynomials may be viewed as the discrete version 
of Hermite polynomials (see, e.g., pj). 

• Statistics. Among the statistics literature we note particularly Ea- 
gleson [To] and Vere- Jones [25] . 



5(3) 



1 1 

3 1 
3 -1 
1 -1 



1 1 

-1 -3 

-1 3 

1 -1 



10 

3 

3 

1 
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• Combinatorics and coding theory. Krawtchouk polynomials are 
essential in Mac Williams' theorem on weight enumerators [121 EI] , and 
are a fundamental example in association schemes [5|, [6l [7] . 

• Probability theory. In the context of the classical symmetric ran- 
dom walk, it is recognized that Krawtchouk's polynomials are elemen- 
tary symmetric functions in variables taking values ±1. It turns out 
that the generating function (11.11) is a martingale in the parameter 

m- 

• Quantum theory. Krawtchouk matrices interpreted as operators 
give rise to two new interpretations in the context of both classical and 
quantum random walks [12]. The significance of the latter interpreta- 
tion lies at the basis of quantum computing. 

Let us proceed to show the relationship between Krawtchouk matrices 
and Sylvester-Hadamard matrices. 

2 Krawtchouk matrices from Hadamard ma- 
trices 

Taking the Kronecker (tensor) product of the initial matrix 



with itself N times defines the family of Sylvester-Hadamard matrices. 
(For a review of Hadamard matrices, see Yarlagadda and Hershey [27].) 

Notation 2.1 Denote the Sylvester-Hadamard matrices, tensor (Kronecker) 
powers of the fundamental matrix H , by 

N times 
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The first three Sylvester- Hadamard matrices are: 



• O • O 

• •00 

• CO* 



• o«o«o«o 

• •©©••oo 

• oo^^oo^ 

• •••oooo 

• o«oo«o« 

• •oooo«« 

• oo«o««o 



where, to emphasize the patterns, we use • for 1 and o for -1. See Table 3 of 
the Appendix for these matrices up to order 5. 

For N = 1, the Hadamard matrix coincides with the Krawtchouk matrix: 
^(1) ^ Now we wish to see how the two classes of matrices are re- 

lated for higher A^. It turns out that appropriately contracting (condensing) 
Hadamard-Sylvester matrices yields corresponding symmetric Krawtchouk 
matrices. 

The problem is that the tensor products disperse the columns and rows 
that have to be summed up to do the contraction. We need to identify the 
right sets of indices. 



Definition 2.2 Define the binary sh.ufB.ing function as the function 

w: N ^ N 

giving the "binary weight" of an integer. That is, let n = dk2^ be the 
binary expansion of the number n. Then w{n) = d^, the number of ones 
in the representation. 

Notice that, as sets, 

w{{0,l,...X-l})-{0,l,...,N} 

Here are the first 16 values of w listed for the integers running from through 
2^ - 1 = 15: 

0112122312232334 
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The shuffling function can be defined recursively. Set 
w{0) = and 

w{2^ + k) = w{k) + 1 (2.1) 

for < k < 2^ . One can thus create the sequence of values of the shuffling 
function by starting with and then appending to the current string of values 
a copy of itself with values increased by 1: 

^ 01 ^ 0112 ^ 01121223 ^ ... 

Now we can state the result; 

Theorem 2.3 Symmetric Krawtchouk matrices are reductions of Hadamard 
matrices as follows: 

sip = E 

w(b)=j 



Example. Let us see the transformation for H^^^ — > 5"^^^ (recall that • 
stands for 1, and o for —1). Applying the binary shuffling function to H^^\ 
mark the rows and columns accordingly: 



0112122312232334 
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The contraction is performed by summing columns with the same index, 
then summing rows in similar fashion. One checks from the given matrix 
that indeed this procedure gives the symmetric Krawtchouk matrix S^^^: 








1 


2 


3 


4 





/I 


4 


6 


4 


1\ 


1 


4 


8 





-8 


-4 


2 


6 





-12 





6 


3 


4 


-8 





8 


-4 


4 




-4 


6 


-4 





Now we give a method for transforming the A^*^ (symmetric) Krawtchouk 
matrix into the N + 1***. 

Definition 2.4 The square contraction r(M) of a 2n x 2n matrix M^b, 
1 < a, 6 < 2n, is the (n + 1) x (n + 1) matrix with entries 



(rM) 



a — 2i, 2i+l 
6=2j,' 2j + l 



ab 



< i,j < n, where the values of M^b with a or 6 outside of the range 
(1, . . . , 2n) are taken as zero. 

Theorem 2.5 Symmetric Krawtchouk matrices satisfy: 
with S^^^ = H. 

Example. Start with symmetric Krawtchouk matrix of order 2: 

"1 2 1 

^(2) =20-2 

1-2 1 

Take the tensor product with H: 

"112 2 11 
1-1 2-2 1-1 

2 2 -2 -2 

2 -2 -2 2 
1 1-2-2 1 1 

1-1-2 2 1-1 
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surround with zeros and contract: 
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Corollary 2.6 Krawtchouk matrices satisfy: 

^(iv+i) _ ^(i^W^W ^ //)(5(iv+i))-i 

where B is the diagonal binomial matrix. 

Note that starting with the 2x2 identity matrix, /, set I^^^ = I, 

Next, we present the algebraic structure underlying these remarkable 
properties. 

3 Krawtchouk matrices and symmetric ten- 
sors 

Given a d-dimensional vector space V over R, one may construct a d^- 
dimensional space V^^ , the A'"-fold tensor product of V, and, as well, a 
('^"'"j^"^) -dimensional symmetric tensor space T/®"^. There is a natural map 

symm: V^^ — > V^^^ 

which, for homogeneous tensors, is defined via 

symm {v <S> w <S> ■ ■ ■) — symmetrization oi {v <S> w <^ . . .) 

For computational purposes, it is convenient to use the fact that the sym- 
metric tensor space of order of a rf- dimensional vector space is isomorphic 
to the space of polynomials in d variables homogeneous of degree N. 
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Let {ei, 62, . . . Cd} be a basis of V. Map Cj to Xi, replace tensor products 
by multiplication of the variables, and extend by linearity For example, 

2ei <8) 62 + 3e2 <8) ei - Tea ® 62 — 5x1X2 - 7x2X3 

thus identifying basis (elementary) tensors in V'^^ that are equivalent under 
any permutation. 

This map induces a map on certain hnear operators. Suppose A e End(y ) 
is a linear transformation on V. This induces a linear transformation Aj^ — 
A'^^ e End(l/®^) defined on elementary tensors by: 

An{v (8)^(8)...)= A{v) (g) A{w) (g) . . . 

Similarly, a linear operator on the symmetric tensor spaces is induced so that 
the following diagram commutes: 



symm 



symm 



This can be understood by examining the action on polynomials. We call A^ 

the symmetric representation of A in degree N. Denote the matrix elements 
of An by A^nn- If A has matrix entries Aij, let 



Vi — Aij X-' 
j 

It is convenient to label variables with indices from to 5 = d — 1. Then the 
matrix elements of the symmetric representation are defined by the expan- 
sion: 

7/"^o . . . n'^s — \^ 'A • • • 

yo ys ~ / J -^mn Xq X^ 



■'S 



with multi-indices m and n homogeneous of degree N. 

Mapping to the symmetric representation is an algebra homomorphism, 

i-e., 

AB^AB 
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Explicitly, in matrix notation, (AB)^^ = ^ (A)mr (-B)rn • 

r 

Now we are ready to state our result 

Proposition 3.1 For each N > 0, the symmetric representation of the N^^ 
Sylvester- Hadamard matrix equals the transposed N^^ Krawtchouk matrix: 



N)ij - J^ji 



Proof: Writing {x,y) for {xo,xi), we have in degree for the k^^ 
component: 

(x + yf-\x-y)' = J2H, 



Substituting x = 1 yields the generating function (11. ip for the Krawtchouk 
matrices with the coefficient of ?/' equal to k\^\ Thus the result. □ 



Insight into these correspondences can be gained by splitting the fun- 
damental Hadamard matrix H (= K^^^) into two special symmetric 2x2 
operators: 

" 1 



1 

1 



G 



so that 



H = F 
One can readily check that 



G 




1 

-1 



F2 = G2 
FH = HG and 



GH = HF 



(3.1) 



The first of the second pair of equations may be viewed as the spectral de- 
composition of F and we can interpret the Hadamard matrix as diagonalizing 
F into G. Taking transposes gives the second equation of (13. ip . 
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Now we proceed to the interpretation leading to a symmetric Bernoulli 
quantum random walk ([I2])- For this interpretation, the Hilbert space of 
states is represented by the N^^ tensor power of the original 2-dimensional 
space V, that is, by the 2^-dimensional Hilbert space . Define the 

following linear operator on \/®^: 

+ . . . 

+/ O / ® ■ ■ ■ O F 

= fl + f2 + ... + fi + ... + fN 

each term describing a "flip" at the i*^ position (cf. [El [22]). Analogously, 
we define: 

Xg = G®/®---®/ 
+ ... 

+/ ® / ® ■ • • ® G 

= 9i+ 92 + ■■■ + 9i + ■■■ + 9n 

From equations (13.1 1) we see that our X-operators intertwine the Sylvester- 
Hadamard matrices: 

= H^^'^Xg and XgH^""^ = H^^'^Xp 

Since products are preserved in the process of passing to the symmetric tensor 
space, we get 

XpHjy = HjsjXg and XgH^ = Hj^Xp (3-2) 

the bars indicating the corresponding induced maps. 

We have seen in Proposition 13. ll how to calculate iJjv from the action of H 
on polynomials in degree N . For symmetric tensors we have the components 
in degree namely x^~^y^, for < < X, where for convenience we 
write X for xq and y for xi. Now consider the generating function for the 
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elementary symmetric functions in the quantum variables fj. This is the 
A'^-fold tensor power 

jr^(t) = (/ + tF)®^ = 7®^ + + • • • 

noting that the coefficient of t is Xp- Similarly, define 

g^{t) ^ (J + tGf^ = 7®^ + + • • • 

From (7 + tF)H = 77(7 + tG) we have 

:r^77(^) = 77(^) and Tn Hn = Hn Qn 

The difficulty is to calculate the action on the symmetric tensors for opera- 
tors, such as Xi?, that are not pure tensor powers. However, from ^iv(t) and 
^Ar(t) we can recover Xp and Xq via 



(7 + tF) 



Xg=- 



t=0 



dt 



(7 + to) 



t=0 



with corresponding relations for the barred operators. Calculating on poly- 
nomials yields the desired results as follows. 



I + tF 



1 t 
t 1 



I + tG 



1 + t 
1-t 



In degree N, using x and y as variables, we get the A;*^'^ component for Xp 
and Xq via 



d 
~dt 



{x + ty)''-\tx + yf = (TV - A;) + k x^-^'^-^^y^"^ 



and since 7 -|- is diagonal. 



d_ 
dt 



(i + i)^-'=(i-t)'=x^-y 



{N - 2k) x^'Y 



t=o 
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For example, calculations for = 4 result in 



X, 
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2 





-2 


-1 


1 





-2 
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-2 
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-1 


1 


-4 


6 


-4 


1 



(3.3) 



(3.4) 



(3.5) 



Since Xq is the result of diagonalizing Xp^ we observe that 



Corollary 3.2 The spectrum of Xp is N, N — 2, ... ,2 — N, —N, coinciding 
with the support of the classical random walk. 



Remark on the shuffling map. Notice that the top row of (/ + tF)®^ 
is exactly t'^^''\ where w{k) is the binary shuffling function of section ^ 
Each time one tensors with I + tF, the original top row is reproduced, then 
concatenated with a replica of itself modified in that each entry picks up a 
factor of t (compare with equation fl2.ll) ). And, collapsing to the symmetric 
tensor space, the top row will have entries {^)t''. This follows as well by 
direct calculation of the 0*^ component matrix elements in degree A^, namely 
by expanding {x + ty)^. 

We continue with some areas where Krawtchouk polynomials/matrices 
play a role, very often not explicitly recognized in the original contexts. 
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4 Ehrenfest urn model 



In order to explain how the apparent irreversibihty of the second law of 
thermodynamics arises from reversible statistical physics, the Ehrenfests in- 
troduced a so-called urn model, variations of which have been considered by 
many authors ([IHl [16], EE]). 

We have an urn with N balls. Each ball can be in two states represented 
by, say, being lead or gold. At each time /c e N, a ball is drawn at random, 
changed by a Midas-like touch into the opposite state (gold lead) and 
placed back in the urn. The question is of course about the distribution of 
states — and this leads to Krawtchouk matrices. 



Represent the states of the model by vectors in M""*"^, namely by the state 
of k gold balls by 



Vfc = [ ■ ■ ■ 1 ■ ■ • 

T 

th 

k position 
In the case of, say, = 3, we have 4 states 



gold balls 
3 lead balls 



1 






1 gold ball 

2 lead balls 




1 





(4.1) 



3 gold balls 
lead balls 






1 



It is easy to see that the matrix of elementary state change in this case is 



1 

3 



2 
3 







1 







^ 

3 ^ 

1 
1 

3 







1 

3 



1 

3 

2 







2 

3 

1 



and in general, we have the Kac matrix with off-diagonals in arithmetic 
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progression 1,2,3,... descending and ascending, respectively: 





N 



1 





2 



- 1 
N -2 





















3 


















■•. 



1 



It turns out that the spectral properties of the Kac matrix involve Krawtchouk 
matrices, namely, the collective solution to the eigenvalue problem Av = Xv 
is 

where A^^^ is the (A^+1) x (A^ + 1) diagonal matrix with entries A^f'' = N — 2i 



N 



N -2 



(*) 



A^-4 



N 



-N 



the (*)'s denoting blocks of zeros. 



To illustrate, for A^ = 3 we have 
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To see this in general, we note that, cf. equations fl3.3H3.5p . these are the 
same operators appearing in the quantum random walk model, namely, we 
discover that A^^) = Xg, A'^^^ = Xj. Now, recalling K^^^ = hJ^, taking 
transposes in equation (13.21) yields 

A(^) ir(^) = ir(^) A(^) and i^(^) A(^) = A(^) iT^^) 
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which is the spectral analysis of A^^^ from both the left and the right. Thus, 
e.g., the columns of the Krawtchouk matrix are eigenvectors of the Ehrenfest 
model with balls where the column := (K.k) has corresponding 
eigenvalue = (A^ - 2k) /N. 

Remarks 

1. Clearly, the Ehrenfest urn problem can be expressed in other terms. For 
instance, it can be reformulated as a random walk on an A^- dimensional 
cube. Suppose an ant walks on the cube, choosing at random an edge 
to progress to the next vertex. Represent the states by vectors in 
Z = Z2 X ■ ■ ■ X Z2, factors. The equivalence of the two problems 
comes via the correspondence of states 

Z 3[aia2...aN] — > G M^+^ 
where w = is the weight of the vector calculated in N, see (14.11) . 

2. The urn model in the appropriate limit as A^ — * 00 leads to a diffusion 
model on the line, the discrete distributions converging to the diffusion 
densities. See Kac' article ([T5]). 

3. There is a rather unexpected connection of the urn model with finite- 
dimensional representations of the Lie algebra s/(2) = so(2, 1). Indeed, 
introduce a new matrix by the commutator: 

A=l[A,A] 

The matrix A is a skew-symmetric version of A. For A^ = 3, it is 

"0-1 " 
^_ 3 0-2 
2 -3 
10 

It turns out that the triple A, A and A is closed under commutation, 
thus forms a Lie algebra, namely 

span {A,A,A} = so{2, 1) ^ s/(2,M) 
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with commutation relations 

[A, A] = 2A , [A, A] = 2yl , [A, A] = -2A 



5 Krawtchouk matrices and classical random 
walks 

In this section we wiU give a probabiUstic meaning to the Krawtchouk ma- 
trices and illustrate some connections with classical random walks. 

5.1 Bernoulli random walk 

Let Xi be independent symmetric Bernoulli random variables taking values 
±1. Let Xn = Xi + ■ ■ ■ + Xn be the associated random walk starting from 
0. Now observe that the generating function of the elementary symmetric 
functions in the Xi is a martingale, in fact a discrete exponential martingale: 

N 

i=l k 

where denotes the A;*'^ elementary symmetric function. The martingale 
property is immediate since each X^ has mean 0. Refining the notation by 
setting a^^^ to denote the k^^ elementary symmetric function in the variables 
Xi, . . . , Xn, multiplying by 1 + vX^^i yields the recurrence 

which, with the boundary conditions af^ = 0, for A; > 0, Cq"^ = 1 for all 
n > 0, yields, for A; > 0, 

N 

(JV+1) _ Y^„(i) y 

that is, these are discrete or prototypical iterated stochastic integrals and 
thus the simplest example of Wiener's homogeneous chaoses. 
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Suppose that at time A^, the number of the Xj that are equal to —1 is j^, 
with the rest equal to +1. Then j^- = {N — XAr)/2 and Mjy can be expressed 
solely in terms of N and xn, or, equivalently, of N and Jn 

= (1 + vf-'^'il - VY'' = (1 + v)iN+-N)/2^^ _ ^)(iV-x^)/2 

From the generating function for the Krawtchouk matrices, equation 
ffO) . follows 

i 

SO that as functions on the Bernoulli space, each sequence of random vari- 
ables K^j^ is a martingale. 

Now we can derive two basic recurrences. From a given column of K^^^ , 
to get the corresponding column in K^^~^^\ we have the Pascal's triangle 
recurrence: 

This follows in the probabilistic setting by writing M^r+i = (1 + vX^)Mn 
and remarking that for j to remain constant, X^v must take the value +1. 
The martingale property is more interesting in the present context. We have 

^ijN - ^(^iiiv+J^i' • • • ,X^) - - (^K^j.^+i + K^j^ j 

since half the time Xat+i is —1, increasing Jn by 1, and half the time Jn is 
unchanged. Thus, writing j for j^, 

which may be considered as a 'reverse Pascal'. 
5.1.1 Orthogonality 

As noted above — here with a slightly simplified notation — it is natural to 
use variables (x, X), with x denoting the position of the random walk after X 
steps. Writing Ka{x,N) for the Krawtchouk polynomials in these variables, 
cf. equation (11.21) . we have the generating function 

TV 

G{v) = ^t;"X,(x,X) = (1 + _ 

a=0 
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The expansion 

il-v)y-%l-{l-R)v)-y = J2lTi<^)n2FA~y r) (5.1) 

n=0 ^ ' 

with (a)„ = r(a + n)/r(a), yields the identification as hypergeometric func- 
tions 

, .r^ fN\ ^ f-aJx-N)/2 

The calculation 

{G{v) G{w)) = JJ(1 + {v + w)Xj + vwX]) = (1 + vwf 

exhibits the orthogonahty of the Ka if one observes that after taking expec- 
tations only terms in the product vw remain. Thus, the are notable for 
two important features: 

1. They are the iterated integrals (sums) of the Bernoulli process. 

2. They are orthogonal polynomials with respect to the binomial distri- 
bution. 

5.2 Multivariate Krawtchouk polynomials 

The probabilistic approach may be carried out for general finite probability 
spaces. Fix an integer d > and d values {C,o, ■ ■ ■ ,^s}, with the convention 
5 — d — 1. Take a sequence of independent identically distributed random 
variables having distribution P(X = ^j) = pj, < j < S. Denote the mean 
and variance of the Xi by /i and as usual. 

For A'^ > 0, we have the martingale 

N 

Mj, = l[{l + v{Xj-fi)) 
We now switch to the multiplicities as variables. Set 

N 

k=l 
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the number of times the value is taken. Thus the generating function 



6 N 



G{v) = J](l + v{^, - = J2 ^"^aK, ...,ns) 

j=0 a=0 

defines our generahzed Krawtchouk polynomials. One quickly gets 

Proposition 5.1 Denoting the multi-index n = (no, . . . ,ns) and by ej the 
standard basis on TL^ , Krawtchouk polynomials satisfy the recurrence 

K^{n + ej) = K^{n) + (^j - /i)if„_i(n) 

We also find by binomial expansion 

Proposition 5.2 

|k|=a j ^ 

where |k| = 'Y^kj. 

j=0 

There is an interesting connection with the multivariate hypergeomet- 
ric functions of Appell and Lauricella. The Lauricella polynomials Fb are 
defined by 



B 



-r, b 

t 



-r)k(b) k 
W|k|k! 



keN« 



with, e.g., r = (ri,...,r5), (r)k = {ri)kr{r2)k2 - ■ ■ {r5)ks for multi-index k, 
also s*^ = s^^ ■ ■ ■ s^*, and k! = /ci! ■ ■ ■ ks\ . Note that t is a single variable. The 
generating function of interest here is 

j reN* ^ ^ 

(5.2) 

a multivariate version of (15. ip . 

Proposition 5.3 Let = |n|. If^Q = 0, then, 



/-r, -n 


1 


■•-) 


V -N 


Pi'' 
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Proof Let Vj = vpj$,j, bj = —rij, t = —N, Sj = Pj in fl5.2p . for I < j < S. 

Note that ^ Vj = vfi, ~ — ^ ~ i Yl '^j) — ''^o- □ 

i<i<<5 

Orthogonality follows similar to the binomial case: 

Proposition 5.4 The Krawtchouk polynomials Ka{no, . . . , ns) are orthogo- 
nal with respect to the induced multinomial distribution. In fact, 
with = |n|, 

Proof 

{G{v)G{w)) = j2( ^ Vo°---^^'n(i+(^+^)(^^-/^)+^^(^^-/^)') 

= [2_^ (Pj + + w)pj{^j - ^) + vwpji^j - 

Thus, {G{v) G{w)) = (1 + vwa'^)'^ . This shows orthogonality and yields 
the squared norms as well. □ 

6 "Kravchukiana" or the World of Krawtchouk 
Polynomials 

About the year 1995, we held a seminar on Krawtchouk polynomials at 
Southern Illinois University. As we continued, we found more and more 
properties and connections with various areas of mathematics. 

Eventually, by the year 2000 the theory of quantum computing had been 
developing with serious interest in the possibility of implementation, at the 
present time of MUCH interest. Sure enough, right in the middle of ev- 
erything there are our flip operators, su(2), etc., etc. — same ingredients 
making up the Krawtchouk universe. Well, we can only report that how this 
all fits together is still quite open. Of special note is the idea of a hardware 
implementation of a Krawtchouk transform. A beginning in this direction 
may be found in the just-published article with Schott, Botros, and Yang [3]. 
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At any rate, for the present we list below the topics which are central to 
our program. They are the basis of the Krawtchouk Encyclopedia, still 
in development; we are in the process of filling in the blanks. An extensive 
web resource for Krawtchouk polynomials we recommend is Zelenkov's site: 

http : //www . geocities . com/ orthpol/ 

Note that we do not mention work in areas less familiar to us, notably 
that relating to g-Krawtchouk polynomials, such as in [23] . 

We welcome contributions. If you wish either to send a reference to your 
paper (s) on Krawtchouk polynomials or contribute an article, please contact 
one of us ! 

Our email: pfeinsil@math.siu.edu or jkocik@math.siu.edu. 
6.1 Krawtchouk Encyclopedia 

Here is a list of topics currently in the Krawtchouk Encyclopedia. 

1. Pascal's Triangle 

2. Random Walks 

• Path integrals 

• A, K, and A 

• Nonsymmetric Walks 

• Symmetric Krawtchouk matrices and binomial expectations 

3. Urn Model 

• Markov chains 

• Initial and invariant distributions 

4. Symmetric Functions. Energy 
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• Elementary symmetric functions and determinants 

• Traces on Grassman algebras 

5. Martingales 

• Iterated integrals 

• Orthogonal functionals 

• Krawtchouk polynomials and multinomial distribution 

6. Lie algebras and Krawtchouk polynomials 

• so (2,1) explained 

• so (2,1) spinors 

• Quaternions and Clifford algebras 

• S and so(2,l) tensors 

• Three-dimensional simple Lie algebras 

7. Lie Groups. Reflections 

• Reflections 

• Krawtchouk matrices as group elements 

8. Representations 

• Splitting formula 

• Hilbert space structure 

9. Quantum Probability and Tensor Algebra 

• Flip operator and quantum random walk 

• Krawtchouk matrices as eigenvectors 

• Trace formulas. MacMahon's Theorem 

• Chebyshev polynomials 

10. Heisenberg Algebra 
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• Representations of the Heisenberg algebra 

• Raising and velocity operator. Number operator 

• Evolution structure. Hamiltonian. 

• Time-zero polynomials 

11. Central Limit Theorem 

• Hermite polynomials 

• Discrete stochastic differential equations 

12. Clebsch-Gordan Coefficients 

• Clebsch-Gordan coefficients and Krawtchouk polynomials 

• Racah coefficients 

13. Orthogonal Polynomials 

• Three-term recurrence in terms of A, K, Lambda 

• Nonsymmetric case 

14. Krawtchouk Transforms 

• Orthogonal transformation associated to K 

• Exponential function in Krawtchouk basis 

• Krawtchouk transform 

15. Hypergeometric Functions 

• Krawtchouk polynomials as hypergeometric functions 

• Addition formulas 

16. Symmetric Krawtchouk Matrices 

• The matrix T 

• S-squared and trace formulas 

• Spectrum of S 
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17. Gaussian Quadrature 

• Zeros of Krawtchouk polynomials 

• Gaussian-Krawtchouk summation 

18. Coding Theory 

• Mac Williams' theorem 

• Association schemes 

19. Appendices 

• K and S matrices for N from 1 to 14 

• Krawtchouk polynomials in the variables x,N/i,j/j,N for N from 1 
to 20 

• Eigenvalues of S 

• Remarks on the multivariate case 

• Time-zero polynomials 

• Mikhail Philippovitch Krawtchouk: a biographical sketch 
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7.1 Krawtchouk matrices 



= [1] 



1 


1 


1 " 








2 





-2 








1 


-1 


1 








1 


1 


1 


1 " 






3 


1 


-1 


-3 






3 


-1 


-1 


3 






1 


-1 


1 


-1 






1 


1 


1 


1 


1 " 




4 


2 





-2 


-4 




6 





-2 





6 




4 


-2 





2 


-4 




1 


-1 


1 


-1 


1 




1 


1 


1 


1 


1 


1 


5 


3 


1 


-1 


-3 


-5 


10 


2 


-2 


-2 


2 


10 


10 


-2 


-2 


2 


2 


-10 


5 


-3 


1 


1 


-3 


5 


1 


-1 


1 


-1 


1 


-1 


1 


1 


1 


1 


1 


1 


6 


4 


2 





-2 


-4 


15 


5 


-1 


-3 


-1 


5 


20 





-4 





4 





15 


-5 


-1 


3 


-1 


-5 


6 


-4 


2 





-2 


4 


1 


-1 


1 


-1 


1 


-1 



1 

-6 
15 

-20 

15 
-6 
1 
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7.2 Symmetric Krawtchouk matrices 



= [1] 
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-1 
i 








i 


— i 








1 


2 


1 " 






2 





-2 






1 


-2 


1 






1 


3 


3 


1 " 




3 


3 


-3 


-3 




3 


-3 


-3 


3 




1 


-3 


3 


-1 




1 


4 


6 


4 


1 


4 


8 





-8 


-4 


6 





-12 





6 


4 


-8 





8 


-4 


1 


-4 


6 


-4 


1 



1 


5 


10 


10 


5 


1 




5 


15 


10 


-10 


-15 


-5 




10 


10 


-20 


-20 


10 


10 




10 


-10 


-20 


20 


10 


-10 




5 


-15 


10 


10 


-15 


5 




1 


-5 


10 


-10 


5 


-1 




1 


6 


15 


20 


15 


6 


1 


6 


24 


30 





-30 


-24 


-6 


15 


30 


-15 


-60 


-15 


30 


15 


20 





-60 





60 





-20 


15 


-30 


-15 


60 


-15 


-30 


15 


6 


-24 


30 





-30 


24 


-6 


1 


-6 


15 


-20 


15 


-6 


1 



Table 2 



Sylvester-Hadamard matrices 



= [•] 
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o 


o 


• 


• 


o 


• 


o 


o 


• 



Table 3 



Replace • with 1 and o with —1 to obtain Sylvester-Hadamard matrices. 



References 



[1] N.M. Atakishiyev and K.B. Wolf, Fractional Fourier- Kravchuk trans- 
form J. Opt. Soc. Amer. A, 14 7 (1997) 1467-1477. 

[2] N. Bose, Digital filters: theory and applications, North-Holland, 1985. 

[3] N. Botros, J. Yang, P. Feinsilver, and R. Schott, Hardware Realization 
of Krawtchouk Transform using VHDL Modeling and FPGAs, IEEE 
Transactions on Industrial Electronics, 49 6 (2002)1306-1312. 

[4] W.Y.C. Chen and J.D. Louck, The combinatorics of a class of represen- 
tation functions, Adv. in Math., 140 (1998) 207-236. 

[5] P. Delsarte, Bounds for restricted codes, by linear programming. Philips 
Res. Reports, 27 (1972) 272-289. 

[6] P. Delsarte, Four fundamental parameters of a code and their combina- 
torial significance. Info. & Control, 23 (1973) 407-438. 

[7] P. Delsarte, An algebraic approach to the association schemes of coding 
theory, Philips Research Reports Supplements, No. 10, 1973. 

[8] C.F. Dunkl, A Krawtchouk polynomial addition theorem and wreath 
products of symmetric groups, Indiana Univ. Math. J., 25 (1976) 335- 
358. 

[9] C.F. Dunkl and D.F. Ramirez, Krawtchouk polynomials and the sym- 
metrization of hypergraphs, SIAM J. Math. Anal., 5 (1974) 351-366. 

[10] G.K. Eagelson, A characterization theorem for positive definite se- 
quences of the Krawtchouk polynomials, Australian J. Stat, 11 (1969) 
29-38. 

[11] P. Feinsilver and R. Fitzgerald, The spectrum of symmetric Krawtchouk 
matrices, Lin. Alg. & Appl., 235 (1996) 121-139. 

[12] P. Feinsilver and J. Kocik, Krawtchouk matrices from classical and quan- 
tum random walks. Contemporary Mathematics, 287 (2001) 83-96. 

[13] P. Feinsilver and R. Schott, Krawtchouk polynomials and fimte probabil- 
ity theory. Probability Measures on Groups X, Plenum (1991) 129-135. 



[14] F.G. Hess, Alternative solution to the Ehrenfest problem, Amer. Math. 
Monthly, 61 (1954) 323-328. 

[15] M. Kac, Random Walks and the theory of Brownian motion, Amer. 
Math. Monthly 54, 369-391, 1947. 

[16] S. Karlin, J. McGregor, Ehrenfest Urn Model, J. Appl. Prob. 2, 352-376, 
1965. 

[17] M. Krawtchouk, Sur une generalisation des polynomes d'Hermite, 
Comptes Rendus, 189 (1929) 620-622. 

[18] M. Krawtchouk, Sur la distribution des racines des polynomes orthogo- 
naux, Comptes Rendus, 196 (1933) 739-741. 

[19] V.I. Levenstein, Krawtchouk polynomials and universal bounds for codes 
and design in Hamming spaces, IEEE Transactions on Information The- 
ory, 41 5 (1995) 1303-1321. 

[20] S.J. Lomonaco, Jr., A Rosetta Stone for quantum mechanics with an 
introduction to quantum computation, 
[http: / / www.arXiv.org/ abs/quant-ph/0007045, 

[21] F.J. MacWilliams and N.J. A. Sloane, The theory of Error- Correcting 
Codes, The Netherlands, North Holland, 1977. 

[22] A.J.F. Siegert, On the approach to statistical equilibrium, Phys. Rev., 
76 (1949), 1708-1714. 

[23] D. Stanton, Some q-Krawtchouk polynomials on Chevalley groups, 
Amer. J. Math., 102 (1980) 625-662. 

[24] G. Szego, Orthogonal Polynomials, Colloquium Publications, Vol. 23, 
New York, AMS, revised eddition 1959, 35-37. 

[25] D. Vere- Jones, Finite bivariate distributions and semi-groups of nonneg- 
ative matrices, Q. J. Math. Oxford, 22 2 (1971) 247-270. 

[26] M. Voit, Asymptotic distributions for the Ehrenfest urn and related ran- 
dom walks, J. Appl. Probab., 33 (1996) 340-356. 



[27] R.K. Rao Yarlagadda and J.E. Hershey, Hadamard matrix analysis and 
synthesis: with applications to communications and signal/image pro- 
cessing, Kluwer Academic Publishers, 1997. 



